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THE  FREE  ELECTRON  LASER  SIDEBAND  INSTABILITY  RECONSIDERED 


I.  INTRODUCTION  AND  SUMMARY 


The  growth  of  parasitic  modes  at  frequencies  near  the  main 


signal  frequency  during  high  power  FEL  operation  was  theoretically 


predicted1’  in  the  early  1980' s.  Since  then  there  has  been  ample 
3  4  5  6 

numerical0’  and  experimental0’  evidence  of  sideband  excitation  in 


constant  wiggler  FELs.  Unstable  modes  in  variable  wiggler  FELs  have 


also  been  observed  in  simulations^-^  and  recently  in  experiment^. 


Sidebands  degrade  the  main  signal  efficiency  and  optical  quality  by 


channeling  a  considerable  fraction  of  the  power  into  parasitic 


frequencies.  The  performance  of  the  mirrors  in  an  oscillator  can  be 


harmed  from  the  modulation  of  the  wave  envelope  caused  by  the 


sidebands.  Last,  but  not  least,  interaction  among  nearby  sidebands 


above  a  certain  amplitude  may  lead  to  chaotic  particle  motion,  loss  of 


trapping  and  incoherent  radiation. 


The  above  have  stimulated  a  considerable  amount  of  theoretical 


work  focused  on  sideband  growth.  Simple  one-dimensional  configurations 


that  are  analytically  tractable  have  been  used  to  model  the  situation. 


Two  lines  of  approach  have  been  considered.  The  single  particle 


picture  regards  the  particle  trajectories  as  functions  of  the  initial 


conditions  and  computes  the  growth  by  ensemble  averaging  over  initial 


distributions  .  The  alternative  approach  assumes  some  adiabatic 


equilibrium  between  the  particles  and  the  main  signal  and  examines  the 


stability  of  the  perturbations  induced  by  the  sidebands,  solving  the 


11  12 

kinetic  equation  ’  .  Because  of  the  equilibrium  assumption  the 
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kinetic  method  is  more  appropriate  for  FELs  operating  at  low  gain. 

In  both  treatments  so  far,  analytic  results  have  been  obtained 
only  for  particles  localized  near  the  bottom  of  the  ponderomotive 
well.  This  implies  the  following  limitations:  The  sideband  spectrum 
becomes  discrete 

«s  =  “r  4  <kr/kw)nV0)’  kr/kw=  2V  (1) 

where  co^(0)  is  the  bounce  frequency  at  the  bottom  of  the  ponderomotive 

well,  k^  k^  are  the  radiation  and  wiggler  wave  numbers  respectively 
2  2  -1/2 

and  y  =(l-v  /c  )  .  The  contribution  from  untrapped  particles  and 
trapped  particles  away  from  the  bottom  is  neglected.  The  effect  of  the 
shear  dw^/dJ,  where  the  action  J  parametrizes  the  distance  from  the 
centre  of  the  separatrix  J=0,  is  lost.  Finally  no  predictions  can  be 
made  about  the  saturation  level  of  potentially  unstable  modes. 

Here  canonical  formalism  is  introduced  by  expressing  the 
unperturbed  particle  orbits  in  terms  of  action-angle  variables.  The 
unperturbed  orbits  are  the  fast  time  averaged  "synchrotron" 
oscillations  of  the  electrons  in  the  potential  well  formed  by  the 
combined  action  of  the  wiggler  and  the  radiation  signal.  The  perturbed 
kinetic  equation  is  solved  in  action  space,  starting  from  an 
equilibrium  extending  over  all  trapped  and  untrapped  electrons.  It  is 
found  that: 

(a)  the  spectrum  becomes  continuous  replacing  c*^( 0 )  by  (^(J) 
in  Eq.  (1).  The  modes  located  at  the  peaks  of  the  unstable  spectrum 
grow  faster,  emerging  as  the  discrete  spectrum  that  is  observed  in 


simulations . 


(b)  More  than  one  groups  of  particles  are  in  resonance  with  a 

*■ 

given  sideband  frequency  «s,  through  different  harmonics  of  their 
bounce  frequency,  and  contribute  to  the  growth  rate. 

(c)  Upper  and  lower  sidebands  located  symmetrically  around  the 
main  signal  frequency  have  opposite  growths  (complementary  stability). 
Therefore  one  mode  is  always  unstable.  There  is  no  stable  distribution 
fQ(J)  except  the  trivial  one  dfQ/dJ=0. 

(d)  The  shear  dw^/dJ  is  stabilizing.  Distributions  with 
gradients  dfQ/dJ  localized  near  the  separatrix  are  found  to  have  the 
minimum  growth  rates  because  of  the  high  shear  there.  This  type  of 
distribution  is  relevant  to  FELs  with  tapered  wigglers. 

(e)  The  growth  is  proportional  to  [df (J)/dw^( J) J ,  the  relative 
population  in  oscillation  quanta  around  resonance,  in  agreement  with 
the  quantum  mechanical  interpretation. 

(f)  For  any  smooth  distribution,  of  finite  dfQ/dJ,  electrons 
at  the  bottom  of  the  well  have  a  negligible  effect  on  stability. 

(g)  Previous  results,  finding  lower  sidebands  having  an 
inherently  larger  growth  than  upper  sidebands,  are  relevant  only  to 
the  limiting  case  of  a  singular  6-function  distribution  fQ(J)=6(J). 
This  case  is  unrealistic  because  a  wide,  smooth  initial  distribution 
in  action  f q( J )  corresponds  to  even  an  ideal  cold  beam  distribution  in 
momentum  fQ(p)=6(p-po) . 

The  nonlinear  saturation  levels  for  the  unstable  modes  and  the 
amplitude  for  stochastic  transition  can  also  be  derived  from  our 
formalism  and  will  be  addressed  in  future  work. 


II.  COMPUTATION  OP  THE  GROWTH 


We  consider  a  circularly  polarized  wiggler  of  constant 
wavelength  2it/kw  and  constant  amplitude  Ay  and  circular  monochromatic 
waves  for  the  main  signal  and  the  sideband.  The  total  vector  potential 
is  given  by 


■  IK-1 


ik  z  i(k  z-w  t) 

ley)Awe  W  -(ex+iey)Are  C  C  -(ex+iey)Ase 


i<ksz-V)j 


CC 


(2) 


where  the  subscripts  w,  r  and  s  stand  for  wiggler,  main  radiation  and 

sideband  respectively,  A^,  Ag,  are  constant  and  Ar>>Ag.  Dispersive 

2 

effects  are  generally  very  small,  of  order  (Wp/w^.)  <<  1,  and  all 

electromagnetic  waves  propagate  with  the  speed  of  light.  Corrections 
to  the  real  part  of  the  wave  number  k(«s>  =  (Wc  are  therefore 

negligible  while  the  imaginary  part  of  kg  can  be  computed  from  the 
wave-particle  energy  balance  equation  without  resorting  to  the 
solution  of  the  dispersion  relation.  Non  linear  phase  slippage  effects 
are  also  omitted. 

The  motion  of  the  electrons  in  the  fields  of  Eq.  (2)  is 
described  by  the  Hamiltonian 


k  -o  ava  a  a 

H(r,<l»;z)  =  —  y  -  — —(cos*  -  *  sin*  )  -  scos(*  -  5  z)  (3) 
rr  Yr  r  rr  s 

where  y  =  y  -  Yr,  time  t  is  normalized  to  of*,  length  z  to  kf mass 

2 

to  m  and  a.  =  eA./m  c  .  Although  k  =  w  =  1  in  the  above  units  we 
e  lie  r  r 

write  them  down  explicitely  to  avoid  some  confusion.  The  resonant 

222  1/2 

energy  y  is  given  by  y  =(l+a  +a  +P  )  with  P  =y  0  and  the 

r  r  v  r  z  z  r  z 
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synchronous  velocity  g  =k  /(k  +k  ).  It  follows  that  the  relation 

z  it  r  v 


between  y_  and  y  is 
r  z 


/i  2* -1/2 

Yz  *  Yr(1  +  aw  > 


Hamiltonian  (3)  was  obtained  for  small  excursions  y/yr  <<  1 
which  is  true  for  electrons  not  too  far  outside  the  separatrix.  It  is 
also  assumed  that  the  wiggler  parameters  change  slowly  compared  to  the 
wiggler  wavelength  2n/kw.  The  term  sin\j»r  parametrizes  the  rate  of 
change  for  the  resonant  energy  caused  by  the  change  in  the  wiggler 
wavelength 


k  a  a 
r  w  r 


where  <J>r=0  corresponds  to  an  untapered  wiggler.  The  terms  proportional 

to  a  a  and  a  a  are  the  ponderomotive  potentials  due  to  the  combined 
w  r  vs 

action  of  the  wiggler  with  the  main  signal  and  the  sideband 
respectively.  The  term  &s  in  Eq.  (4)  signifies  the  Doppler  shifted 
frequency  departure  of  the  sideband  from  the  main  signal 


*s  =  <kw  /kr>  (ks"  V* 


When  as=0  the  Hamiltonian  Hq(P,^)  is  exactly  integrable.  The 
unperturbed  trajectories,  shown  in  Fig.  1(a),  are  given  by 


where  the  constant  K,  regarded  as  the  reduced  energy  in  the 
ponderomotive  frame,  is  determined  from  the  initial  conditions 


K=Hq(Po,  .  These  trajectories  take  the  simplest  possible  form 
expressed  in  terms  of  the  action  angle  variables  (J,6)  defined  as 


J  =  ^  <fd*  y (K,*),  ©  =  fjj^d^  y(K,*')  ,  (8) 

where  K  =  Hq(J)  and  the  path  of  integration  is  over  the  unperturbed 
orbits.  Clearly  the  action  J  is  related  to  the  area  in  phase  space 
enclosed  by  the  orbit.  For  untrapped  particles  and  in  case  of 
untapered  wiggler  the  limits  of  <|/  integration  are  from  0  to  2ji.  In 
case  of  a  tapered  wiggler  we  may  define  the  limits  of  integration 
starting  at  <J/s  =  2n  -  \J»r  and  going  back  to  \l*s  around  the  separatrix. 
Thus  J  remains  finite,  avoiding  an  infinite  jump  in  action  across  the 
separatrix  that  would  otherwise  result  by  considering  the  full  orbit 
length  for  unbound  orbits.  In  the  above  variables  the  width  of  the 
separatriv  is  given  by  Jg,  the  value  of  action  at  the  separatrix. 

We  now  apply  the  transformation  (8)  to  the  Hamiltonian  Eq.  (3) 
resulting  into  the  decomposition  of  the  perturbation  cos[  <H  J , 0)- &gt ] 
into  harmonics  of  the  synchrotron  oscillation  and  obtaining 


H(J,9)  =  H  (J)  +  ?n{  h+(J)cos(n0+S  t)  +  h  (J)cos(n0-&  t)  }, 
o  n — u  n  s  n  s 


h±(j)  = 


a  a 
w  r 


Q*(J). 


(9) 


where  Qn±(J)  are  the  Fourier  coefficients  for  the  synchrotron 
harmonics.  This  expansion  is  similar  in  spirit  with  the  more  familiar 


case  of  expanding  a  plane  wave  propagating  across  a  magnetized  plasma 
into  cyclotron  harmonics.  Expression  (9)  for  the  Hamiltonian  is 
formally  independent  on  the  choice  of  tapered  or  untapered  viggler 
which  will  only  affect  the  functional  forms  for  Qn±(J). 

Here  we  present  the  analysis  for  constant  parameter  wiggler 
were  the  transformations  are  solvable  in  closed  forms.  The  relation 
between  new  (J,8)  and  old  (P,\|i)  variables  is  expressed  by 


Js  [E2(X)-(1-XZ)E1(X)] , 


J  = 


2Jg  X  E2(l/X), 


X2<1 


2 

XZ>1 


,  a  a  k  .1/2 

Js=  HsP 

w 


Sln2=i  ,1 


X  snl^CXie], 
snlifijd/Xie], 


2 

X'Xl 


XZ>1 


(10) 


with  £^*^2  ^e  complete  elliptic  integrals  of  the  first  and  second 
kind  and  sn  the  Jacobi  elliptic  sine  function.  The  trapping  parameter 
X2  is  given  by 


X2,  K  t  G 


a  a 
w  r 


2G  ’ 


(11) 


2  2 

and  we  have  X  <1,  X  >1  for  trapped  and  untrapped  particles 

2 

respectively.  The  three  constants  of  the  motion  X  ,  K  and  J  are 
mutually  dependent  and  any  one  of  them  defines  uniquely  a  trajectory. 

Since  Hq(J)  is  independent  of  9  the  unperturbed  trajectories 
in  J,  9  are  straight  lines 


7 


J(z)  =  const . , 


0(z)  =  eQ  +  Kb(J)z,  (12) 

Kb(J)  =  dHQ(J)/dJ. 

<b(J)  is  the  synchrotron  wavenumber,  connected  to  the  synchrotron 
length  Lb  by  «b  =  2n/Lb  and  to  the  synchrotron  frequency  J )  by 


VJ> 


VJ) 

*zc  ' 


From  (8), (12)  and  (13)  we  find 


(13) 


/  «b(0)2E™(X)'  X  <X 

VJ)  ={  *  2  (14> 

V0)ifn7x>’  x  >x 

where 

c^z  1/2 

V°)  -  ~T  <awarW  * 

r 


The  coefficients  Qn±(J)  are  obtained  by  contour  integration  in  the 
complex  plane  using  the  double  periodicity  properties  of  the  Jacobi 
elliptic  functions, 


(si) 


n  nil 


n 

2 


Ej(X) 


l-(-q) 


X2<1 


(15) 


2,2 

nn  X 


Ej(l/X) 


■"  ( 


itE^l/X), 


1-q 


2n  1 


1+q 


2n 


O'  'V  \ 

’  o  =  exP  l  Ejli/xTJ 


X2<1 


Qn(J)  tends  to  0  as  J  tends  to  either  0  or  ®  while  the  maximum  occurs 


at  the  separatrix  X  =1.  Qn(J)  becomes  progressively  smaller  with 


increasing  n  everywhere  except  in  the  vicinity  of  the  separatricies. 


That  limits  the  number  of  harmonics  to  be  considered  for  instabilities 


localized  away  from  the  separatrix.  The  perturbed  equations  of  motion 


in  the  standard  form  read 


J  =n?on{h^(J)sin(n0+6sz)  +  hn(J)sin(n@  -$sz)}, 


■X  Uil  Uli 

dl°  =  VJ)  +n!o{dTnc0s(ne+Ssz)  +  djncos(n0-6sz)}- 


The  growth  gg  for  a  given  sideband  is  determined  by  the  energy 


balance  equation 


h  9  9  An  i  ik  z-ieo  t 

df  V  2gsas=  “17“  {  <  jy(z>t)  l(iase  +  cc)  >  +  cc  }  (17) 

s  } 


The  bracket  <  >  denotes  the  average  of  a  quantity  <f>  =  L  J  f(z)  dz 


over  the  short  length  scale  L  =  2n/kw  <<  =  2ji/k ^ ,  and  j  is  the 


transverse  current 


jy(z,t)  =  enb 


p2  rt  r® 

d0  dJ  v 

J0  Jq  * 


Sf(0, J;z) 


Vy  =  Y  py  =  -(e/ ymc)Ay(z , t ) , 


t(z)  =  z/c&z- 


For  small  sideband  signal  6f  is  the  perturbed  distribution 


under  the  Hamiltonian  flow,  Eqs.  (16) 


■A.  /.V  V  SJ 


(19) 


!fi  K(*_h  dJdA 

3z  +  VJ'3G  "  '  dz  dJ  ’ 

where  f  is  the  equilibrium  distribution  along  the  unperturbed  orbits 
dfQ/dz=0.  Any  distribution  fQ(J)  depending  on  J  alone  is  invariant 
under  Eqs.  (12) 


,  j  t  d  f 

d_£  _  dJ_  _ o 

dz  o  _  dz  dJ 


=  0 


thus  f  is  not  uniquely  defined. 

Ue  now  assume  a  small  imaginary  part  in  the  sideband 
wavenumber  =  Im(ks).  To  obtain  gg  for  a  given  equilibrium  f q ( J )  we 
solve  Eq.  (19)  for  f.,  then  substitute  Eqs.  (9)  and  (18)  into  (17)  and 
finally  average  over  the  short  length  scales  L.  The  computation  is 
similar  to  that  described  in  Ref.  13,  with  the  exception  that  here  the 
equations  of  motion  and  the  evolution  of  the  distribution  function  are 
parametrized  by  the  length  z  rather  than  the  time  t,  yielding 


-  00  n2oo2  »  |Q+(J  )  I2 

7  2  r  P  r  1  nv 

rdf  i 
0 

r  w  .  2  ",  r 

oo  4y  w  n=l  r 

,dJ  > 

,dJ  > 

s  r  r 

J 

n 

(20) 


with  Wp  the  beam  plasma  frequency.  The  superscript  +  or  -  corresponds 

to  upper  w  >w  and  lower  oo  <oo  sideband  respectively.  J  is  given 
s  r  sr  r  J  n  ° 


implicitly  by  the  resonant  condition 


±nwb(Jn)  =  (kw/kr)  (W’ 


(21) 


10 


as  illustrated  in  Fig.  1(b).  The  sum  over  n  on  the  right-hand  side  of 
Eq.  (20)  includes  the  contribution  from  all  resonant  groups  of 
particles.  The  action  Jn  in  Eq.  (17)  labels  the  orbit  having  the  nth 
harmonic  of  the  local  bounce  frequency  in  resonance  with  the  sideband. 
The  growth  is  determined  by  the  slope  of  the  distribution  function  f 
near  these  resonant  orbits  J=Jn. 

III.  RESULTS  AND  CONCLUSIONS 

The  opposite  signs  in  the  right-hand  side  of  Eq.  (20)  denote 

that  upper  and  lower  sidebands  located  symmetrically  around  the  main 

signal  |(os+-«r|  =  |ws  |  have  opposite  growth  rates  >  g+/g  =  -1  + 
2 

0(6w/«r)  .  The  physical  reason  is  connected  with  quantum  mechanical 

considerations.  The  practical  implication  is  that  some  sideband  is 

always  unstable  except  in  case  of  trivial  equilibrium  dfQ/dJ=0.  It 

also  reveals  that,  once  destabilized,  upper  and  lower  modes  have 

growth  rates  of  equal  magnitude  and  are  equally  dangerous.  We  clarify 

that,  in  general,  the  opposite  signs  do  not  imply  that  all  upper 

sidebands  have  the  same  kind  of  stability,  opposite  to  the  stability 

of  all  lower  sidebands.  Depending  upon  fQ  stability  may  change- sign 

between  two  upper  (or  lower)  frequencies  because  the  slopes  dfo/dJn 

change  as  the  location  of  the  resonant  J  's  shift  with  w  . 

ns 

For  a  monotonic  distribution  f  we  observe  that  trapped  and 
untrapped  particles  yield  opposite  contributions  to  a  given  mode 
because  dw^/dJ  changes  sign  across  the  separatrix.  If  trapped 
particles  are  stabilizing  untrapped  are  destabilizing  and  vice  versa. 


High  shear  is  stabilizing  tending  to  reduce  the  magnitude  of  the 
growth  gg.  This  is  expected  as  the  number  of  resonant  orbits  is 
inversely  proportional  to  |d«^/dJj.  Shear  tends  to  infinity  near  the 
separatrix  thus  the  modes  generated  by  electrons  orbiting  near  the 
separatrix,  corresponding  to  small  lws“wrl  <<  (k^/k^nto^O) ,  have  the 
smallest  growth  rates.  Modes  coming  from  electrons  near  the  center 
Jn=0  with  |ws-«r|  =  (kr/kv)nw^(0)  have  also  small  growth  rates  as  the 
coefficients  Qn(J)  tend  to  zero  there.  In  particular,  given  any  smooth 
distribution  f q ( J ) ,  electrons  exactly  at  the  bottom  of  the  well  have  a 
null  contribution  to  the  instability.  It  takes  a  singular  distribution 
of  infinite  gradient  at  J=0  to  create  an  instability  from  electrons 
exactly  at  the  center  of  the  separatrix.  The  shear  is  generally  higher 
for  untrapped  particles;  thus  the  trapped  ones  usually  dominate  the 
instability. 

Although  Eq.  (20)  was  obtained  on  purely  classical  arguments 

it  nevertheless  admits  the  correct  quantum  mechanical  interpretation. 

Expressing  (dfQ/dJ)(d<»^/dJ)_^  as  dfQ/dw^  we  observe  that  the  growth  gs 

is  proportional  to  the  difference  in  population  of  oscillation  quanta 

between  the  energy  levels  across  the  resonance.  A  more  detailed 

explanation  for  the  opposite  signs  in  Eq.  (20)  is  given  in  Ref.  13. 

The  normalized  growth  is  plotted  against  the  percentage 

mismatch  (ws-wj.)/«i_  for  both  upper  and  lower  sidebands  in  Figs.  2  and 

4.  The  contribution  up  to  the  third  harmonic  n<4  in  Eq .  (20)  is 

included  in  these  plots.  The  parameters  chosen  correspond  to  a  wiggler 

_4 

wavelength  A  =3  cm,  a  =5,  main  signal  strength  a  =5x10  ,  beam  energy 

w  w  s 

2 

of  16.50  MeV  (y  =  32.3)  and  current  density  j=100  A/cm  (beam  density 
10  -3 

6.25x10  cm  ).  The  equlibrium  distribution  is  a  Gaussian 
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?  1/9  9  2 

f0(J)*»(l/2nD  exp(-(J-JQ)V2D  )  centered  halfway  inside  the  island, 

Jq=Js/2,  and  of  width  D  equal  to  half  the  separatrix  distance  D>Jg/2. 

We  plot  the  contribution  of  only  the  fundamental,  n=l  in  Eq.  (20),  in 

Fig.  2(a)  adding  the  first  harmonic  n=2  in  (b)  the  second  harmonic  in 

(c)  and  the  third  harmonic  in  (d).  New  unstable  bands  emerge  with  each 

harmonic,  while  the  growth  for  already  unstable  bands  is  modified.  For 

example  we  observe  two  upper  and  two  lower  unstable  bands  in  (b)  but 

only  one  upper  and  three  lower  bands  in  (d).  We  find  the  contributions 

from  higher  than  the  third  harmonics  n>4  generally  negligible.  The 

upper  frequency  ws>“r  and  the  lower  frequency  ws<wr  parts  of  the 

unstable  spectrum  come  from  the  regimes  of  negative  and  positive 

slope  dfQ/dWk  respectively,  shown  in  Fig.  3.  The  lower  sideband  growth 

is  peaking  at  frequencies  corresponding  to  nwjj(Jmax)  with  Jmax  the 

value  maximizing  d^/dw^J).  The  peaks  for  the  upper  sideband  growth, 

however,  do  not  occur  at  J=0  that  minimizes  df  /dco^J),  but  at  J 

halfway  inside  the  negative  slope  regime.  This  is  because  Qn(J)  and 

consequently  g  are  zero  at  J=0,  showing  the  negligible  contribution 
s  . 

from  electrons  at  the  bottom  of  the  well.  The  analogous  effect  in 
plasma  physics  is  the  elimination  of  the  thermal  effects  when  goes 
to  zero.  In  any  case  the  most  unstable  modes  are  far  from  the 
frequencies  |  =  (kr/kv)n«^(0)  pointed  by  the  arrows  in  Fig.  2(d). 

Since  sidebands  can  not  be  shut  out  completely  it  remains 
debatable  whether  a  distribution  function  can  be  tailored 
experimentally  to  minimize  their  growth  rate.  From  the  previous 
discussion  a  flat  distribution  inside  the  trapped  regime  with  sharp 
gradients  localized  at  the  separatrix  seems  the  appropriate  choice. 
Instabilities  will  then  localize  near  the  separatrix  and  the  growth 
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will  be  suppressed  by  the  strong  shear.  To  check  this  we  plot  the 

growth  rates  from  Eq.  (20)  for  two  types  of  distributions  f Q ( J ) :  (i) 

2  1/2  2  2 

Two  Gaussians  fQ( J)=(l/2nD  )  exp(-J  /2D  )  centered  at  the  centre  of 

the  island  and  of  characteristic  lengths  D  equal  to  half  the  island 

width  D=J  /2  in  Fig.  4(a)  and  one  island  width  D=J  in  Fig.  4(b).  (ii) 

s  s 

Two  "step-like"  distributions  of  the  form  f0(J)*(l/aD)*^exp(-(J/aD)N) 

1/N 

with  N=16.  Selecting  a»(N/N-l)  places  the  sharp  gradient  at  J=D  and 

we  plot  the  case  D=J  /2  in  Fig.  4(c)  and  D=J  in  Fig.  4(d). 

s  s 

Comparing  Fig.  4(a)  to  Fig.  4(b)  and  Fig.  4(c)  to  Fig.  4(d)  it 
is  seen  that  the  growth  rates  between  similar  types  of  distributions 
tend  to  decrease  as  the  location  of  the  maximum  gradient  (dfQ/dJ)max 
approaches  the  separatrix.  In  both  cases  there  is  more  than  one  order 
of  magnitude  reduction  in  the  growth  by  shifting  the  maximum  gradient 
position  from  D=J  /2  to  D=J  .  Because  f  was  chosen  monotonic  in  all 
above  plots  and  because  it  was  limited  to  trapped  particles,  df^dc^ 
preserves  sign  and  only  lower  modes  are  unstable.  The  spectral  width 
of  the  unstable  regimes  is  reduced  with  a  parallel  increase  in  the 
maximum  growth  as  one  goes  from  the  Gaussian  type  to  the  step-like 
type  of  distributions.  Also  the  distance  of  the  sideband  frequencies 
from  the  main  signal  decreases  by  shifting  the  gradient  position  D 
closer  to  the  separatrix.  Distributions  with  sharp  gradients  at  the 
separatrix  such  as  those  in  Figs.  4(b)  and  4(d)  are  perhaps  more 
relevant  to  the  case  of  variable  wiggler  FEL,  where  the  "bucket"  of 
the  trapped  particles  is  decelerating  in  phase  space  leaving  the 
untrapped  particles  behind,  and  only  small  diffusion  occurs  across  the 
separatrix  allowing  sharp  gradients. 

The  limit  of  a  6-function  distribution  fo=6(J-Jo),  examined 
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elsewhere  ,  is  the  most  unstable  case  but  of  the  least  practical 
interest,  because  even  the  case  of  a  monoenergetic  beam  distribution 
Pz  =pQ  is  described  in  J-space  by  a  smooth  distribution  fQ(J)  of 
finite  width  AJ  (see  Fig.  1). 

Our  calculations  done  for  the  case  of  a  fixed  wavelength 
wiggler  can  be  easily  extended  to  the  case  of  a  variable  wiggler  FEL 
provided  that  the  same  adiabatic  assumptions  hold.  Only  the  functional 
relations  Eqs.  (7)  between  the  action-angle  variables  (J,0)  and  the 
coordinates  (P,<J>)  need  to  change.  The  derivation  of  the  growth  rate, 
performed  in  J-space,  is  independent  of  the  transformation  J(P,^)  and 
the  result,  Eq.  (20),  stands  as  is.  It  is  the  high  growth  of  tapered 
wiggler  FELs  that  seems  to  challenge  the  adiabatic  approach  to  the 
problem.  In  this  case  the  change  of  the  signal  amplitude  ar(t)  in  time 
together  with  the  dependence  of  the  equilibrium  f  on  both  time  and  0 
should  be  included  for  a  more  realistic  treatment. 
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Figure  1.  Time  averaged  motion  without  the  sidebands.  (a)Plots  in 

phase  space  of  the  unperturbed  orbits  Ho(P,*)-K.  The 

intersections  with  the  horizontal  line  P»const.  mark,  the 

initial  conditions  for  each  orbit,  (b)  The  normalized  bounce 

frequency  and  the  first  two  harmonics  as  functions  of  the 

2 

trapping  parameter  X  (J).  The  intersections  with  the 

horizontal  line  w»(«  -w  )k  /k  determine  the  position  J  of 

s  r  w  r  r  n 

the  resonant  orbits  for  a  given  w  . 


9/wr  Q/cu 


2  2 

Figure  2.  Growth  for  a  Gaussian  distribution  fo(J)«Cexp(-(J-JQ)  /2D  ) 
centered  halfway  inside  the  trapped  particle  island  Jq»=Js/2 
and  of  width  D  equal  to  Jg.  The  normalized  growth  gs/<*>r  is 
plotted  versus  Sot/u^  for  (a)  the  fundamental  contribution  n=l 
in  Eq.  (20),  (b)  including  the  first  harmonic  n-2,  (c)  two 
harmonics  and  (d)  three  harmonics.  The  most  unstable  modes  do 
not  correspond  to  harmonics  of  the  bounce  frequency  at  the 
bottom  <^(0)  indicated  by  the  arrows. 
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Figure  3.  Plot  of  dfo(J)/d«j)(J)=df0/dJ(d(Ojj/dJ) as  a  function  of 
K(J)-K(0)  for  the  distribution  fQ  of  Figure  2.  The  slope 
goes  to  zero  near  the  separatrix  K(Jg)  because  of  the 
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